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Abstract. The mean width is a measure on n-dimensional convex 
bodies. An integral formula for the mean width of a regular n-simplex appeared 
in the electrical engineering literature in 1997. As a consequence, expressions 
for the expected range of a sample of n + 1 normally distributed variables, for 
n < 6, carry over to widths of regular n-simplices. As another consequence, 
precise asymptotics for the mean width become available as n — )• oo. 

Let C be a convex body in M". A width is the distance between a pair of parallel 
C-supporting planes (linear varieties of dimension n ~ 1). Every unit vector m G 
determines a unique such pair of planes orthogonal to u and hence a width w{u). 
Let u be uniformly distributed on the unit sphere S*""^ C M". Then w is a random 
variable and 
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E(w3) = — arccos -- 
for C = the regular 3-simplex (tetrahedron) in M.^ with edges of unit length and 
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for C = the regular 4-simplex in with edges of unit length [ll!2]. Our contribution 
is to extend the preceding mean width results to regular ra-simplices in M" for n < 6. 
We similarly extend the following mean square width result: 

E (wl) = 

which, as far as is known, first appeared in [3]. 

The key observation underlying our work is due to Sun [3], which in turn draws 
upon material in [5l |6] . It does not seem to have been acknowledged in the mathe- 
matics literature. After most of this paper was written, we found [7], which assigns 
priority to to Hadwiger [8] and to Ruben [9J for closely related ideas. 
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1. Order Statistics 
Let Xi, X2, Xn denote a random sample from a Normal (0, 1) distribution, that 
is, with density function / and cumulative distribution F: 

—00 

The first two moments of the range 

Tn = max{Xi,X2, ...,X„} - min{Xi,X2, 
are given by [TOl [TT] 



00 

/^n = E(r„)= y {1-F(x)'^-[1-F(xr}cia;, 



00 y 

= E(r. 



l) = 2j J {1- F{yr - [1 - F{x)r + [F{y) - F{x)r} dx dy. 



6.156..., 



For small n, exact expressions are possible [12 

/i2 = ^ = 1.128..., z/2 = 2, 

/i3 = ^ = 1.692..., ^3 = 2 (1 + ^) = 3.653..., 

/i4 = ^ (1 - 2^2) = 2.058..., u, = 2(l + ^) = 5.012..., 

/.5 = ^ (1 - 3^2) = 2.325..., Z.5 = 2 f 1 + ^ + f - 

/ie = ^ (1 - 4^2 + 2T,) = 2.534..., z/^ = 2 M + '-^^^ - fS^ - ^^S^ = 7.142..., 

^^^(l_5^2 + 5r2) = 2.704..., 7 y , ^ 1/2 ,2 ,3 
' ^ +3|v^T3 + ^f/-^y) = 8.007... 

where 

7r/4 

^ Vk f dx 1 o , IN 

6fc = / — = = — arcsec [k + 1) , 

- ' y/k + sec(x)2 27r ^ ^' 



TT 



7r/4 7r/4 tt 

^ /• /• 1 / , fc(fc + l) , ^ 

y + sec(x)2 + sec(y)2 " 27r2 7 + - tan(z)^ ' 

00^ 
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1 1 
jj _ I f arcsec {2t^ + 4) , t/ _ ^ /" arcsec {t^ + 5) 

(2^2 + 1) v/2t2T3 7727 {t^ j^2) V'FTl 



7r2 



The preceding table complements an analogous table in [15] for first and second 
moments of max{Xi, X2, Similar expressions for /ig = 2.847... and = 
8.778... remain to be found. 

2. Key Observation 
Let us rescale length so that the circumradius of the ra-simplex is 1. Adjusted width 
will be denoted by Wn- Using optimality properties of the n-simplex, Sun |4j deduced 
a formula for mean half width: 



n + 1 /(n + l)n T (f ) f^f^y-^.,, 

-E(wn) = \- tHtt F(^] f(x)dx 

2 ^ ' 2 V 27r r (^^) J ^ ' 



[n 



+ ip r(f) 



2n r(^) 



xF{xyf{x) dx 



(see Corollary 2 on p. 1581 and its proof on p. 1585; his M is the same as our n + 1). 
We recognize the latter integral as fi^^i/ (2(n + 1)); hence 



+ 1 r (-) 
^K) = \l 2n r(4^) ^"+^ 

and therefore 

1 r (-) 
EK) = -^rp|fy/^n+i 

because, in our original scaling, the circumradius is \Jn/{2{n + 1)). 

No similar integral expression for E (w^) appears in [4J. We circumvent this diffi- 
culty by noticing that the formula [Ml Hj 




72. 



r(f) 



bears some resemblance to the coefficient of in our expression for E (w„). The 
square version 

/ n 

-2 



n 
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is trivial and leads us to conjecture that 

by analogy. Numerical confirmation for n < 6 is possible via the computer algebra 
technique described in [3]. 

In summary, we have mean width results 

E{w2) = - = 0.954929658551372..., 

71 

E (ws) = ^ (1 - 2S2) = 0.912260171954089..., 
20 

E{w4) = — (1 - 3^2) = 0.874843256085440..., 

Stt 

45 

E (W5) = — (1 - 4^2 + 2T2) = 0.842274297659162..., 
16 

56 

E (we) = — (1 - 5^2 + 5T2) = 0.813743951590337..., 

57r 

and mean square width results 

E (wl) = ^ I 1 + ^ 1 = 0.913496671566344..., 
2 \ zn I 

E (wl) = \(^^ ^"V^) " 0.835419517991054..., 
w( ^\ 1 A^5v^ 30 5v^ 

4 \ ZIT 71 71 

= 0.769572883591771..., 

. 2x 1 , 5(9 + 2v/3) 90 15v^ 

= 0.714241915072694..., 
.2. 1 / 35v^ 210 105 35v^ 

+ ^T3 + ^f/-l^V 
= 0.667314714095430.... 
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3. ASYMPTOTICS 

We turn now to the asymptotic distribution of r„ as n — )■ oo. Define a„ to be the 
positive solution of the equation [121 IE] 



2vr al exp (a^) 

that is, 



271 J " ' ' 2 v^21n(n) 

in terms of the Lambert W function [19]. It can be proved that the required density 
is a convolution 



lim — P ( a/2 Infn) (r„ — 2a„) < y] = / exp(— x — e~^) exp(— — x) — e~^^~''^)rfx 
n^oo ay \ /J 

— oo 

= 2 e-^i^o (2 e-^/') 

where Kq is the modified Bessel function of the second kind [22] . A random variable 
Y, distributed as such, satisfies 

E(r) = 27, E(F2) = y + 

where 7 is the Euler-Mascheroni constant [23j. This implies that 

J 7 ^ o ATTTT ln(lnH)+ln(47r)-27 

/i„ ~ 2 a„ + , ^ ^ ~ 2V21n(n) 



v/2 \n{n) J v^2 ln(n) 



and hence 



2n V 4ny V 2n ln(n) J y ^/2h^{ 
2 /ln(n) ln(ln(n)) + ln(47r) - 27 



n 



n 2A/nln(n) 
More terms in the asymptotic expansion are possible. 
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If we rescale length so that the inradius of the n-simplex is 1 and denote adjusted 
width by Wn, then 



E(wn) ~ V2n ■ 2\ ~ 2J2n\n(n) 

V n 

because, in our original scaling, the inradius is A/l/(2n(n + 1)). This first-order 
approximation is consistent with [2]. 

4. Regular Octahedron 
As an aside, we return to the setting of and review our computational methods 
for C = the regular octahedron with edges of unit length. 
For simplicity, let () be the octahedron with vertices 

t;i = (1,0,0), ^72 = (-1,0,0), t;3 = (0,1,0), 
V4 = (0, -1, 0) , v, = (0, 0, 1) , ve = (0, 0, -1) . 

At the end, it will be necessary to normalize by -\/2, the edge-length of 0- 

Also let be the union of six overlapping balls of radius 1/2 centered at fi/2, 
V2/2, V3/2, V4/2, Vq/2, Vq/2. Clearly C and has centroid (0,0,0). A diameter 
of is the length of the intersection between and a line passing through the origin. 

Computing all widths of is equivalent to computing all diameters of (>• The 
latter is achieved as follows. Fix a point (a, b, c) on the unit sphere. The line L 
passing through (0, 0, 0) and (a, b, c) has parametric representation 

X = ta, y = tb, z = tc, tGM 

and hence y = {b/a)x, z = {c/a)x assuming a 7^ 0. The nontrivial intersection 
between first sphere and L satisfies 

thus Xi = since a'^ + b'^ + = 1; the nontrivial intersection between second sphere 
and L satisfies 

(x+i)V(^x)V(fx)^ = l 

thus X2 = —a^. The nontrivial intersection between third/fourth sphere and L 
satisfies 

thus = a6 , 0:4 = —ab. The nontrivial intersection between fifth/sixth sphere and 
L satisfies 
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thus X5 — a c , xe — —a c. 

We now examine all pairwise distances, squared, between the six intersection 
points: 



(X, 



= < 







2 






if i = i,j = 


2 






1 - 2a 6 - 


if i = = 


3 or i = 


2,i 


= 4 


1 + 2a 6 - 


if i = L j = 


4 or i = 




= 3 


1 - 2a c - 6^ 


if i = l^j = 


5 or i = 


2,J 


= 6 


1 + 2a c - 52 


if i = i,j = 


6 or i = 


2,J 


= 5 




if i = 3, j = 


4 






{b - cf 


if i = 3, j = 


5 or i = 


4,J 


= 6 


{b + cf 


if i = 3, j = 


6 or i = 




= 5 


if i = 5, j = 


6 







and define 

g{a,b) 



max {4a^ 1 - 2a 6 - c^, 1 + 2a 6 - c^, 1 - 2a c - 6^ 1 + 2a c - 6^ 
4bMb-cf,{b + c)\4c^}. 



The mean width for C is 

27r 



— =— / / \/ g{cos 9 sin (/?, sin 6 sin (/?, cos </?) sin (/? = — arccos ( - ) 
V2 47rj J TT \3/ 





and the mean square width is 



27r 7r 



-— I / 5'(cos 6* sin sin 9 sin cos smipdipd9 = - (1 

2 47ry J 3 \ 

00 ^ 

Here are details on the final integral. A plot of the surface 



2^3' 



TT 



g{cos 9 sin 93, sin 9 sin cp, cos (f) 



appears in Figure 1, where < 6' < 27r and < 93 < tt. Figure 2 contains the same 
surface, but viewed from above. Our focus will be on the part of the surface to the 
right of the bottom center, specifically < 9 < n /4 and 7r/2<(/7^9/4. The volume 
under this part is 1/24*^ of the volume under the full surface. 



Mean Width of a Regular Simplex 



8 



We need to find the precise upper bound on as a function of 9. Recall the formula 
for g as a maximum over nine terms; let ge denote the i^^ term, where 1 < i < 9. 
Then the upper bound on is found by solving the equation 

gi (cos 9 simp, sin 9 simp, cos ip) = gg (cos 6* sin yj, sin 9 simp, cos (p) 

for ip. We obtain ip{9) = 2 a.Tctan{h{9)) , where 



/ 3 + cos(2g) 
h{9) = cosb* + Y — - 

and, in particular, 

ip{0) = 2 arctan (l + ^ 2.3562, 

(^(7r/4) = 2 arctan (^(^1 + Vtj /V2^ ^ 2.1862. 
It follows that g = gi foT < 9 < TT /A and tt/2 < ip < 2 arctan(/i). Now we have 



gi (cos 9 sin ip, sin ^ sin ip, cos ip) sin ip dip 



1 1 
24^, 

(1 + cos(2^)) (cos(3(/?) - 9 cos(¥?)) 



and 



487r 



x|2arctan(h) (1 + 4/l + /l^) (1 - 4/l + /l^) (1 - /l2^ 



therefore 



/2 (1 + /i2)3 

/ X |2arctan(h) 1 — /i^ 

-j^ i'2 arctan(/i) 



2 47r 



(7i (cos 9 sin yj, sin ^ sin yj, cos if) sin c?*/? 

n/2 



{h^ + 4/^2 + 1) (/i2 - 1) (1 + cos(20)) 
67r(l + /i2)3 ■ 

Integrating this expression from to 7r/4 gives the desired formula for E ("U^octa)- 

5. n-CUBES 

After having written the preceding, we discovered [7j, which gives the mean width 
for a regular ra-simplex in M" as 



Mean Width of a Regular Simplex 



9 



Consistency is readily established; nothing is said in f7] about the connection between 
E (wn) and order statistics from a normal distribution (more precisely, the expected 
range 

By contrast, the mean width for an ra-cube with edges of unit length is elementary: 



and we conjecture that 



E (^Icube) = 1 + 



TT 

6. n-CROSSPOLYTOPES 

A regular n-crosspolytope with edges of unit length has mean width [TJ El] 



oo 



2y/2n{n-l) T (f ) . ^ 

l^^n-crosspolytope j ~ -p / n+1 \ I ^ eril^Xj (XX, 

V 2 y g 

the case n = 3 corresponds to the octahedron discussed earlier. It is not surprising 
that a connection exists with order statistics from a half-normal (folded) distribution. 
We will examine this later, as well as relevant expressions from An appropriate 
mean square conjecture also needs to be formulated in this scenario. 
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Figure 1: Surface plot of \fgj^-, where < ^ < 27r and < </? < tt. 



Figure 2: Another view of y/gj^, with contours of intersection. 



This figure "Figure01.jpg" is available in "jpg" format from: 
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This figure "Figure02.jpg" is available in "jpg" format from: 
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